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C/> ABSTRACT 

Context. The nearby Type la SN 201 lfe has provided an unprecedented opportunity to derive some of the properties of the progenitor system, 
I which is one of the key open problems in the supernova (SN) field. 

Aims. This study attempts to establish a link between the reasonably well known nature of the progenitor and its surrounding environment. This 
is done with the aim of enabling the identification of similar systems in the vast majority of the cases, when distance and epoch of discovery do 
not allow a direct approach. 

O^l Methods. To study the circumstellar environment of SN 201 lfe we have obtained high-resolution spectroscopy of SN 201 lfe on 12 epochs, 
from 8 to 86 days after the estimated date of explosion, targeting in particular at the time evolution of Can and Nai. 

Results. Three main absorption systems are identified from Can and Nai, one associated to the Milky Way, one probably arising within a 
H high-velocity cloud, and one most likely associated to the halo of M101 . The Galactic and host galaxy reddening, deduced from the integrated 
y—i equivalent widths (EW) of the Nai lines are £ fl _ v =0.01 1±0.002 and £Vv=0.014±0.002 mag, respectively. The host galaxy absorption is 
r- H dominated by a component detected at the same velocity measured from the 21-cm Hi line at the projected SN position (~ 1 80 kms -1 ). 

During the ~3 months covered by our observations, its EW changed by 15.6+6.5 mA. This small variation is shown to be compatible with the 
. geometric effects produced by the rapid SN photosphere expansion coupled to the patchy fractal structure of the ISM. The observed behavior 
is fully consistent with ISM properties similar to those derived for our own Galaxy, with evidences for structures on scales <100 AU. 
Conclusions. SN 201 lfe appears to be surrounded by a "clean" environment. The lack of blue-shifted, time-variant absorption features is fully 
^ consistent with the progenitor being a binary system with a main-sequence, or even another degenerate star. 

Key words, supernovae: general - supernovae: individual: SN 201 lfe - ISM: dust, extinction 



1. Introduction 

Since the introduction of the accreting White Dwarf (WD) sce- 
nario (Wheelan & Iben 1973 ), the question mark about the na- 
ture of the progenitor systems of Type la Supernovae (hereafter 

* Based on observations collected at the Mercator telescope, 
Telescopio Nazionale Galileo, Nordic Optical Telescope at Roque de 
los Muchachos, La Palma (Spain), and at the 1.82 m Copernico tele- 
scope on Mt. Ekar (Asiago, Italy). 



SN la) has been growing bigger and bigger (see Patat l20lT1 foi' a 
general review on the subject). One of the most beaten tracks to 
obtain information on the progenitor system (and more specif- 
ically on the nature of the donor star) has been the search for 
imprints of circumstellar material (CSM) lost by the system 
prior to the explosion. Notwithstanding the extensive attempts 
undertaken to detect such material, as of today no signs of it 
have been found in the form of optical (Lundqvist 2003 2005 ; 
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Table 1. Log of the observations 



Date 


UT 


JD 


Phase 


Epoch 


Airmass 


Exposure 


Instrument 


Hel. Correction 




(start) 


(start) 


(days) w 


(days)'**' 


(start) 


(seconds) 




(kms- 1 ) 


31-08-2011 


20:39 


2455805.357 


-11 


8.2 


1.64 


1800 


FIES 


-4.19 


03-09-2011 


20:20 


2455808.344 


-8 


11.2 


1.52 


3600 


SARG 


-3.46 


07-09-2011 


20:20 


2455812.340 


-4 


15.1 


1.58 


1800 


HERMES 


-2.47 


13-09-2011 


19:23 


2455818.308 


+2 


21.1 


1.60 


4200 


AES 


-0.94 


13-09-2011 


20:16 


2455818.342 


+2 


21.1 


1.72 


1350 


HERMES 


-0.96 


14-09-2011 


18:57 


2455819.290 


+3 


22.1 


1.50 


3600 


AES 


-0.69 


21-09-2011 


18:59 


2455826.291 


+ 10 


29.1 


1.63 


5400 


AES 


+ 1.09 


22-09-2011 


18:40 


2455827.278 


+11 


30.1 


1.56 


9000 


AES 


+ 1.35 


29-09-2011 


20:47 


2455834.363 


+ 18 


37.2 


2.63 


2700 


FIES 


+3.10 


18-10-2011 


17:54 


2455853.246 


+37 


56.1 


1.84 


7200 


AES 


+7.68 


15-11-2011 


16:52 


2455881.203 


+65 


84.0 


2.80 


4800 


AES 


+ 12.86 


18-11-2011 


01:28 


2455883.561 


+67 


86.4 


2.65 


12600 


AES 


+ 13.52 


18-11-2011 


05:05 


2455883.710 


+67 


86.5 


2.92 


4800 


SARG 


+ 13.55 



(*) Phases from V maximum were computed assuming a typical rise time of 19 days. 

(**) Epochs are computed from the estimated explosion time (JD=2455797.196; Nugent et al. 1201 let . 



Mattila et al. 12005b . radio (Panagia et al. 120061 Chomiuk et al. 
|20TTT i, and UV/X-Ray emission (Immler et al. 120061 . 

Independent information on the binary system can be ob- 
tained studying the effects the explosion has on the compan- 
ion star. So far two possibilities have been explored: the col- 
lision of the SN ejecta with the companion star surface, and 
the stripping of the companion's envelope. While the former is 
expected to re-heat part of the SN ejecta causing a luminosity 
excess in the very early epochs (Kasen 1201 Oi l, the latter is sup- 
posed to entrain some of the companion atmospheric material, 
which should eventually become visible in the form of emis- 
sion lines at late epochs (see e.g. Marietta, Burrows & Fryxell 
2000 and references therein). So far, no evidences of either 
phenomena have been found (Hayden et al. 120101 Bianco et al. 
I20TT1 Leonard [2007b . 

These facts constitute a considerable body of evidence in 
favor of SN la progenitor systems where the companion star is 
either a MS star, or even another degenerate body. However, it 
is only with the very recent study of SN 201 lfe that this was di- 
rectly demonstrated with a high level of confidence. SN 201 lfe 
(PTFllkly) was discovered by the Palomar Transient Factory 
on August 24 UT in the outskirts of the nearby (~6 Mpc), face- 
on, SAB(rs)cd galaxy M101 (NGC 5457), and immediately 
classified as a very early Type la (Nugent et al. 1201 lal 1201 lbl >. 
The SN was not detected (down to mag 20.6) on an image taken 
on August 23, implying the object was caught within one day 
after the explosion. 

Based on pre-explosion HST images, Li et al. (1201 11 1 could 
place for the first time a direct and stringent limit on the com- 
panion's luminosity. This rules out a red-giant donor star, and 
hence a symbiotic system as the progenitor for this particular 
event. The observations presented by Li et al. (2011) exclude 
the He-star channel only partially, while a system with a MS- 
star transferring mass via Roche lobe overflow is still consistent 
with the data, as is a double-degenerate configuration. The lack 
of an early shock (Nugent el al. 2011c) goes along the same 
lines, ruling out companion stars with extended atmospheres. 



Thanks to the very early detection, Nugent et al. (1201 lei 
could put a stringent upper limit to the exploding star radius 
(<0.1 Rq), hence providing for the first time a direct evidence 
that the progenitor is indeed a compact star (see also Brown 
et al. 1201 II) . Although a carbon-burning MS star cannot not be 
excluded, a C-O WD is the most favored candidate. The range 
of allowed objects was further narrowed down by Bloom et 
al. d2011| l, who derived a more stringent upper limit for the 
progenitor radius (0.02 R ), and concluded that the exploding 
star must indeed be a WD (for independent considerations on 
the progenitor's nature, see also Liu et al. 1201 II) . Bloom et al. 
(1201 11 1 also place an upper limit to the radius of the companion 
(0.1 R Q ), which excludes Roche-lobe overflowing giants and 
MS stars. Finally, prompt radio and X-ray observations placed 
the most stringent limit ever on the mass loss rate for a SN la: 
dM/dt<10~ 8 (w/100 km s 1 ) M G yr 1 (Horesh et al. l20TT1 l. 

All these facts give a very consistent view on the progenitor 
system of SN 201 lfe, which is very plausibly identified as a C- 
O WD accreting material from a MS star (Li et al. 120 1 It Brown 
et al. l20m Nugent et al. 1201 lcl Bloom et al. [20TTb . This is 
a major break-through in the field. Having studied an object 
whose progenitor system is reasonably well known, we can 
now ask ourselves a new question. Would we recognize similar 
systems in other SN la, whose distance and discovery time do 
not allow the close scrutiny that was possible for SN 201 lfe? 
The key issue here is that one needs to consider the largest 
possible number of observables, because, for instance, similar 
spectra and light curves may come from different channels. 

In this article we focus on the properties of the circumstel- 
lar environment, whose study has led to promising results in the 
last years (Patat et al. 2007a, Simon et al. 2009, Sternberg et 
al. |201 11 1. For this purpose, and following the previous studies, 
we have obtained multi-epoch, high-resolution, optical spec- 
troscopy of SN 201 lfe covering the first three months of its 
evolution. 

The paper, reporting the results of this campaign, is struc- 
tured as follows. In Sect. [2] we give an account of the observa- 
tions and data reduction. Sect.|3]recaps the detection of narrow 
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Table 2. Instrumental setups 



ID 


Site 


Telescope 


Res. 


Disp. 


Range 








power 


(A/pixel) 


(A) 


AES 


Asiago 


1.8-m 


18 000 


0.126 


3800-7500 


FIES 


La Palma 


2.6-m 


48 000 


0.036 


3650-7350 


HERMES 


La Palma 


1.2-m 


82000 


0.016 


3770-9000 


SARG 


La Palma 


3.5-m 


66000 


0.021 


4600-6150 



absorption features, whereas Sect.[4]discusses the implications 
on the extinction suffered by the SN. In Sect. [5] we illustrate 
the properties of the intervening gas as derived from Voigt pro- 
file fittings, and we study their behavior as a function of time. 
The modest time evolution is then analyzed in Sect. [6] in the 
context of the small-scale structure of the inter-stellar medium. 
The results of this study in connection to the progenitor's na- 
ture are finally discussed in Sect. [7] which also summarizes our 
conclusions. 



2. Observations and data reduction 

We obtained high-resolution spectroscopy of SN 201 lfe on 12 
epochs, ranging from 8 to 86 days after the estimated date 
of explosion (JD=2455797.2; Nugent et al. 1201 lcl >. The ob- 
serving log is presented in Table [T] which includes exposure 
times, phase from maximum light, epoch from explosion, air- 
mass, and heliocentric correction. Phases are only indicative, 
and were computed assuming a typical SNIa rise time of 19 
days. The observations were carried out with four different in- 
strumental setups. These are summarized in Table [2] and de- 
tailed in the next subsections. 

2.1. Asiago Echelle Spectrograph 

Seven spectra of SN 201 lfe were obtained with the REOSC 
Asiago Echelle Spectrograph (hereafter AES) attached to the 
1.82 m telescope on Mt. Ekar (Asiago, Italy). The spectrograph 
is equipped with an an E2V CCD42-40 AIMO, back illumi- 
nated detector (2048x2048 pixel, 13.5 fim in size). A 2" slit 
was used, with a fixed east-west orientation, yielding a resolv- 
ing power R - A/AA ~ 18 000. The wavelength range cov- 
ered by the spectra is 3780-7620 A. The data reduction was 
performed with tasks within the IRAfQ ECHELLE package fol- 
lowing the standard procedures. All images were first over- 
scan corrected and trimmed. A normalized flat-field image was 
then derived from observations of a halogen lamp and used to 
flat field the science exposures. Next, each spectral order was 
traced and the spectra were extracted with the optimal extrac- 
tion algorithm of Home ( 1986), simultaneously subtracting the 
scattered light. The transformation from pixels to wavelength 
was derived from observations of a Thorium- Argon arc-lamp. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, under contract with the National Science 
Foundation. 



2.2. FIES 

Two spectra (Aug 31 and Sep 29, 201 1) were obtained with the 
fiber-fed (1.3") echelle spectrograph (FIESj^J mounted at the 
Nordic Optical Telescope (NOT) on La Palma (Spain). We used 
the medium-resolution fiber, that provides a resolving power 
R~48 000, and continuous spectral coverage from 3700 to 7300 
A in a single exposure. Data were reduced using standard IRAF 
tasks, as outlined in the previous section. 

2.3. HERMES 

Two spectra (Sep 7 and 13, 2011) were obtained with the 
High Efficiency and Resolution Mercator Echelle Spectrograph 
(HERMES; Raskin et al. 120111 . mounted on the Mercator 1 .2 m 
telescope in La Palma. HERMES is a high-resolution fiber- 
fed (2.5") echelle spectrograph offering spectral resolution of 
R -80000, and full spectral coverage from 3770 to 9000 A 
in a single exposure. Data were processed with an automated 
data reduction pipeline that includes bias correction from the 
prescan region on the CCD, flat-field correction via Halogen 
lamps, wavelength calibration through Th-Ar reference spec- 
tra, cosmic clipping, and background modelization. 

2.4. SARG 

Two spectra (Sep 3, and Nov 18, 2011) were obtained with 
the Spettrografo Alta Risoluzione Galileo (SARG; Gratton 
et al. 120011 1. mounted at the 3.5 m Telescopio Nazionale 
Galileo (TNG) in La Palma. SARG is a cross-dispersed echelle 
spectrograph which provides spectral resolution ranging from 
R=29 000 up to R= 164,000. For our observations we used the 
cross-disperser #3 grism, a GG455 order sorting filter, and slit 
#2 (0.8"x 5.3"). This setup yields R~66,000 and a wavelength 
range 4600-6150 A. The data reduction was performed using 
standard IRAF tasks, as outlined in Sect. |2. 1| 

2.5. Removal of telluric features 

Although the spectral region of Can H&K is free of telluric 
absorption features, this is not the case for Nai D, and espe- 
cially Ki /17665 and /17699. Because the observations were 
conducted mostly at relatively high airmasses (see Table [TJ 
these are quite severely disturbed by water bands. 

The removal of telluric absorption features is usually per- 
formed using spectra of bright stars obtained with the same 
instrumental setup and similar airmass as for the scientific tar- 
get. In this paper we adopted a different approach, using an 
atmospheric transmission model (see for instance Stevenson 
|1994t Seifahrt et al. 1201 Oi l. The advantages over a direct ob- 
servation are: a) the absence of intrinsic photospheric features, 
b) the absence of interstellar features (chiefly Na i D lines), c) 
the infinite signal-to-noise, and d) the exact wavelength scale. 
The latter provides an absolute reference to which the observed 
spectra can be corrected via cross-correlation techniques (see 
below). 



http://www.not.iac.es/instruments/fies/ 
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Fig. 2. Narrow absorption components towards SN 201 lfe for Can K (upper panel; day 21) and Nai D2 (lower panel; day 1 1.2). 
Main absorption features are marked by vertical dashed lines. 



In this work we adopted the Line By Line Radiation 
Transfer Model (LBLRTM; Clough et al. l2005T l. This code is 
based on the HITRAN database (Rothman et al. 12009b . It has 
been validated against real spectra from the UV to the sub- 
millimeter, and is widely used for the retrieval of atmospheric 
constituents^] LBLRTM solves the radiative transfer using an 
input atmospheric profile, which contains the height profiles 
for temperature, pressure, and chemical composition. For our 
purposes we have adopted the standard LBLRTM equatorial 
profile. 

After computing the synthetic model in the regions of inter- 
est and for the given airmass, this is convolved with the instru- 
mental profile and the wavelength correction is computed via 
cross-correlation. The absence of intrinsic narrow absorptions 
in the SN spectrum and the weakness of the interstellar fea- 
tures in the Na 1 and K 1 regions (see below) make the correla- 
tion function peak very well defined. The wavelength solutions 
for FIES, HERMES and S ARG spectra were all found to match 
the model wavelength with a maximum deviation of 130 m s _1 , 
that is more than 5 times smaller than the pixel size. Therefore, 
no wavelength correction was applied to these data. In the case 
of the AES spectra, maximum deviations of 6.6 km s (about 



See http : //rtweb . aer . com/lblrtm . html 



one pixel) were found, and the spectra were corrected accord- 
ingly. The final wavelength accuracy is ~0.1 km s for FIES, 
HERMES and SARG data, while this is -0.3 km s _1 for the 
AES spectra. 

An example of telluric correction is presented in Fig.[T]for 
the SARG spectrum. All telluric features are removed at the 
level of a few percent, which is comparable to (or better than) 
the continuum RMS noise of all spectra presented in this paper. 

As a for the sky background subtraction, we note that in 
some spectra the terrestrial Na 1 D emission components (both 
natural and artificial) are not completely removed. Because of 
the wavelength shift of most features, this problem does af- 
fect only a Galactic absorption close to null restframe velocity 
(see Sect. [3]l, whose profile (and thus equivalent width) is most 
likely affected by the incomplete sky subtraction. 



3. Narrow absorption features in SN 201 1fe 

3.1. Can H&K and Na i D 

A number of Ca n H,K and Na i D components are detected at 
different velocities, as shown in Fig. [2] which presents the high- 
est resolution spectra available for the two wavelength ranges 
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we derive from our data (EW (Nai D 2 )=47+l mA, EW(Nai 



5897 5898 5899 5900 5901 5902 
Observed Wavelength (A) 

Fig. 1. Example of telluric features removal for the SARG 
spectrum in the region of Nai D 2 . Upper panel: observed data 
(black line) and best fit LBLRTM model (blue line). Lower 
panel: telluric -corrected spectrum. The dashed vertical line 
marks the position of Nai D 2 (vhel ~180 kms -1 ). The signal- 
to-noise ratio on the continuum is about 100 (the error-bars in- 
dicate the estimated RMS noise). 



D!)=27±l mA; see Sect. 5.1 1 is consistent with that reported 
by Nugent et al. , velocity discrepancy is statistically very sig- 
nificant (~14 km s -1 ), and cannot be explained in terms of the 
quoted errors. Once applying the heliocentric corrections listed 
in Table [T] to the different spectra, the velocities deduced for 
the Na I D components all agree to within a few 0.1 kms . 
The average heliocentric velocity deduced from the best fit pro- 
file of the most intense host galaxy component (see Sect. |3j is 
179.6±0.2 km s -1 . Given the very consistent results we obtain 
from different setups and epochs, we believe our result is cor- 
rect to within the reported errors. 

3.2. Other absorptions 

A careful inspection of the telluric-corrected spectra covering 
the relevant wavelength range shows no trace of the two typi- 
cal K i lines at 7665 and 7699 A, both for the Galaxy and for 
M101 (at the velocities derived for the Can and Nai features. 
See previous section). For the host galaxy the 5-cr detection 
limit derived from the HERMES spectrum taken on day 15 is 
EW(Ki ,17699)<5 mA. The Ki ,17665 is severely affected by a 
telluric absorption, and so the upper limit estimate for this line 
is more problematic. However, the two lines have very simi- 
lar strengths, and so the limit derived above applies to both of 
them. Similar considerations apply to the non-detection of Ti n, 
Cai, CHandCH+. 

Finally, none of the known strong Diffuse Interstellar Bands 
(DIBs; 5780, 5797, 6284, 6379, 6613; Herbig [1995]) are de- 
tected. 



(obtained with HERMES and SARG, respectively; see Tables[T] 
and [2). 

Two main Galactic absorption systems are detected around 
-70 and -20 kms 1 (marked as A and B in Fig. [2]), with 
weaker and diffuse absorptions possibly present at intermedi- 
ate velocities (see the Can profile). A third system is clearly 
detected at about +130 kms -1 (C). The recession velocity of 
M101 is 241±2 kms -1 (de Vaucouleurs et al. 119911) . and there- 
fore this component may be generated within a high-velocity 
cloud (HVC) belonging either to the Milky Way or to the host 
galaxy. The radial velocity derived from 21 -cm line observa- 
tions at the projected SN site is about 180 km s _1 (Bosma, 
Goss & Allen 1981 1 their Fig. 1). This makes the association to 
M101 more plausible, although not certain (see also Sect. [5}. 
Finally, we detect an absorption system with two main com- 
ponents at about 180 and 195 km s (D). These velocities are 
fully consistent with the radial velocity field at the projected 
SN location (Bosma et al. [1981). This strongly suggests that 
the corresponding absorptions arise in M101. All features are 
detected in both H and K lines, while system C is not detected 
in Nai. The component close to 195 km s _1 is clearly detected 
in Nai only in the highest signal-to-noise spectrum (SARG, 
day 11) with a total equivalent width EW=7±2 mA, while it is 
only marginally detected in all other spectra. 

Nugent et al. (1201 lcl l report the detection of one single ab- 
sorption feature (EW=45±9 mA), which they identify as Na I 
D 2 at v=194±l km s _1 . Although the weighted average EW 



4. Reddening 

The galactic Nai D components are very weak (see Fig. [2]), in- 
dicating that SN 201 lfe suffers a very low Galactic reddening. 
The integrated EW(Na 1 D 1 ) derived from our highest signal-to- 
noise spectrum (SARG, day 11) for the Galactic components 
(systems A and B in Fig. |2| is 38±5 mA. Applying the re- 
lation inferred by Munari & Z witter (1997) for Galactic red- 
dening, we get E B _ V =0.01 1+0.002 mag, which is consistent 
with the expected Galactic extinction along the line of sight 
to M101 (£Vv=0.009 mag; Schlegel, Finkbeiner & Davies 
1998). According to the Munari & Zwitter relation, the ex- 
pected EW for the Ki ,17699 is ~4 mA, which is fully con- 
sistent with the lack of K 1 detection in our spectra. 

As shown in the previous section, we detect two features 
most probably associated to the host galaxy, at ~180 and 
~200 km s _1 , with the latter contributing only marginally to 
the total Nai column density. The integrated D 2 EW mea- 
sured on the SARG spectrum (day 11) is 47±2 mA (EW(Nai 
D!)=25±l mA). Applying the Munari & Zwitter relation, this 
turns into £Vi/=0.014±0.002 mag (Ay=0.04 mag), in agree- 
ment with the valued derived by Nugent et al. d2011cl) using 
the same method. 

Independent constraints on the reddening can be placed 
studying the DIBs (see for instance Sollerman et al. 2005 Cox 
& Patat 2008). As it turns out, these are much less stringent 
than those derived from the sodium lines, and we mention them 
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Table 3. Na i and Ca n interstellar cloud model parameters 





Weighted Average 


Day 8 


Day 11 


Day 15 


Day 21 


Day 37 


V 


b 


N 


N 


N 


N 


N 


N 


Nai 


-67.2 
-23.3 
-12.6 
179.6 
193.6 


2.1 
3.0 
3.4 
3.8 
6.4 


10 52 +0 08 
10 74+ 07 

10 55 +o n 

1 1 44+ 002 
10 57 +012 

-0.16 


10 52 +018 

1U,J -0.26 

10.70_2| 3 
10.18_£* 

1 1 43 +003 
10.56!^° 


10.53*"* 
10.80!°;^ 
10 64 +0 06 

lu - -0.06 

11.43!g : °! 
10.46!^ 


10.36!™ 
10.8QUJS 

10 27 +0 - 20 

1 1 43 +0 02 
10 50 +0 16 


10 62 +0 0S 

i 7 3+o.o6 
iu. /o_ 008 

10 66 +0 0S 
iu.uu_ 009 

1 1 44+0.02 

11 -^-0.02 

10 64 +0 11 

lU.Ut _„ j, 


10.50_™ 
10.60!™ 
10.80^ 
11 51 +003 
10 42 +0 26 



Can 



-66.4 


0.7 


11 10 +017 


11 05 +034 


11.15_j»f 


11 07 +018 


-65.7 


32.0 


1 1 90 +007 


12 00 +o os 


11.85_™ 


1 1 89 +0 09 

11,0 -0.10 
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here only for the sake of completeness. Assuming a standard 
relation between DIB strength and Eg-v (e.g. Luna et al. 2008 1 
EW(6613)/£ , ij - 1 /=310 mA), and based on the lack of DIB de- 
tection (see previous section) gives an upper limit on the M101 
reddening of 0.4 mag. A similar assessment for the 5780 DIB 
(EW(5780)/£_._ V = 460 mA) gives E(B-V) < 0.43 mag. Due 
to the intrinsic weakness of other DIBs, the upper limits on 
their EW derived from our spectra give even less stringent con- 
straints. For example, the non-detection of the 6379 DIB yields 
£ B -v <0.65 mag (EW(6379)/£_- V =88 mA; Luna et al. 2008). 



5. Narrow absorption line modeling 

In order to derive interstellar cloud velocities, Doppler widths, 
and column densities we modeled the best resolution and 
signal-to-noise spectra using the vapid routine (Howarth et al. 
2002). Saturation effects were implicitly taken into account fit- 
ting simultaneously lines of the same species (Nai Di and D2; 
Can H and K). The results for the FIES, HERMES and SARG 
(first epoch only) spectra are shown in Figs. [3] and |4] for Nai 
and Ca 11. The Na 1 spectra require 5 components for a good fit, 
whereas 9 components were required for Ca 11. Error estimates 
were derived using Monte Carlo noise replication. The mod- 
els were re-fit to 1000 replicated data-sets, each with Gaussian 
noise of RMS equal to that of the continuum added. The ±68% 
ranges in the resulting data give rigorous l<x errors on each pa- 
rameter. Results are given in Table [3] for each epoch, as well 
as the results of a weighted-average fit to all epochs, in which 
the fitted data points were weighted by the continuum signal- 
to-noise ratio of the respective spectra. The data show no ev- 
idence for statistically significant variations in cloud compo- 
nent velocities or Doppler b parameters between epochs (see 
below), so the v and b values were held fixed at the weighted- 
average fit values and only the jV values were allowed to vary 



for each epoch. Column densities (cm 2 ) are given as base-ten 
logarithms. The units of v and b are km s~' . 

The analysis of Milky Way interstellar gas is beyond 
the scope of this paper. Here we mention that for the most 
marked Galactic feature (-23 kms~'), the N(Na i)/N(Ca 11) ra- 
tio is ~0.50, which is typical of galactic halos (Hobbs 1978; 
Crawford [T992l Bertin et al. 119931 . Sembach & Danks (fT994l > 
measured a mean Galactic Halo N(Nai)/N(Can)~0.5. 

Interestingly, the N(Na i)/N(Ca 11) ratio for the main fea- 
ture in M101 (+180 kms -1 ) is -0.43, indicating that it may 
as well originate in the halo of the host galaxy (see also 
King et al. [19951 for the case of SN 1994D). The really de- 
viant case is the one of the component at ~128 kms -1 , which 
is clearly detected in Can but not in Nai. The upper limit 
to the Nai column density deduced from the high-signal to 
noise SARG spectrum (day 1 1) is 1.5xl0 10 cm 2 , which yields 
N(Nai)/N(Can)<0.04. Such low values are rarely observed, 
and indicate a mild depletion of calcium onto dust grains. A 
similar case was observed for one high velocity component 
along the line of sight to SN 2000cx, where Ca 11 was detected 
while Nai was not (Patat et al.__007b). The corresponding ra- 
tio was N(Nai)/N(Can)<0.1. Finally, King et al. (TT995T > found 
N(Na i)/N(Ca n)=0. 1-0.3 for the high-velocity clouds observed 
along the line of sight to SN 1994D. Low ratios are, in gen- 
eral, typical of high velocity gas (Siluk & Silk 1974). The small 
value we derived for the component at ~ 180 kms -1 is similar 
to those observed in intermediate- and high-velocity clouds by 
Molaro et al. ( 1993), who concluded that such clouds may con- 
tain very little dust, with calcium and sodium mainly in gaseous 
form (see also the properties measured by Cox et al. 2007] to- 
wards LMC sight-lines). For these reasons it is reasonable to 
conclude that absorption system C (see Fig.|2]i arises in a HVC. 
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Fig. 4. Ca ii K spectra for each epoch for which these data were 
obtained (black histograms), with best-fit models overlaid (red 
curves). Spectra have been normalised and offset vertically for 
display. Heliocentric velocity scale is given with respect to the 
Ca ii K line rest wavelength. 



Fig. 3. Na i D spectra for each epoch (black histograms), with 
best-fit models overlaid (red curves). Spectra have been nor- 
malised and offset vertically for display. Heliocentric velocity 
scale is given with respect to the Na D 2 line rest wavelength. 
The spurious emissions observed in some spectra around and 
300 km s _1 are residuals of the sky subtraction (see Sect. [2]). 



5.1. Time evolution 

The derived cloud parameters show some marginal evidence 
for time variability. The only possibly significant variation is 
seen for the host galaxy Na 1 component at ~ 1 80 km s~ 1 . While 
the logarithmic column density derived for this absorption is 
11.44+0.02 on days 8, 11, 15 and 21, on day 37 this grows to 
11.52±0.03, which indicates a variation at the 3-sigma level. 
One may be tempted to attribute this fluctuation to the different 
instrumental setups. However, a direct comparison of the two 
spectra taken on days 8 and 37 (both obtained with FIES and 
the same identical setup), clearly shows that both components 
Di and D2 had consistently deepened by about 16% on the sec- 
ond epoch, while their FWHM remained unchanged. The inte- 



grated EWs directly measured on the two spectra, are 23.6+1.9, 
46.1+2.0 mA, and 30.4+2.4, 53.0+2.5 mA for Di and D 2 on 
the two epochs respectively, implying a variation AEW(Nai 
D2)=6.9+3.2 mA. This strengthens the conclusion that what- 
ever the reason for the change is, this cannot be attributed to a 
variation of instrumental resolutior|f] 

A similar check on the EWs can be done using the two 
SARG observations, which were obtained ~75 days apart (see 
Table[T]i. The high airmass, and the limited visibility window of 
SN 201 lfe on the second epoch, did not allow us to reach a very 
high signal-to-noise ratio (~16 at 5900 A), so that the EW esti- 
mate is rather uncertain. However, the difference is AEW(Na 1 
D2)=13.3+6.4 mA, which confirms the variation at the ~2-rr 
level. 

Although AES spectra have a significantly lower resolution 
(16.7 kms~', see Table [2i, given the simple structure of the 



component at ~ 180 km s~ 



they can be used together with the 



4 Being a fiber-fed spectrograph, FIES is not subject to changes of 
resolution due to seeing variations. Also, a check on the comparison 
lamp emission lines on the two epochs shows that the resolution in- 
deed remained constant. 
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Table 4. Na i D equivalent widths 



Setup 


Epoch 


EW(D 2 ) 


EW(D0 


S/N 


Resolution 




(days) 


(mA) 


(mA) 


(5900A) 


(km s- 1 ) 


(*) 


1.5 


45.0±9.0 


- 


- 


6.0 


FIES#1 


8.2 


46.1±2.0 


23.6+1.9 


60 


6.3 


SARG#1 


11.2 


46.6±1.5 


26.3+1.0 


115 


4.5 


HERMES* 1 


15.1 


44.3±1.8 


27.6+1.8 


50 


3.7 


AES#1 


21.1 


49.2±2.2 


34.2+1.9 


130 


16.7 


HERMES#2 


21.1 


44.3±1.8 


27.6+1.8 


50 


3.7 


AES#2 


22.1 


53.8±4.2 


28.5+4.0 


105 


16.7 


AES#3 


29.1 


51.9+2.7 


31.0+2.2 


190 


16.7 


AES#4 


30.1 


44.5±2.3 


33.7+2.5 


155 


16.7 


FIES#2 


37.2 


53.0±2.5 


30.4+2.4 


55 


6.3 


AES#5 


56.1 


48.7±7.1 


34.4+7.3 


50 


16.7 


AES#6+#7 


85.2 


58.7±3.7 


33.5+3.0 


115 


16.7 


SARG#2 


86.4 


59.9+6.2 


23.7+7.3 


16 


4.5 


Weighted average 


47.1+0.8 


27.4+1.1 







(*)The first epoch is from Nugent et al. ( 120 1 lcl 7 



higher resolution data to study the variability on the whole time 
range spanned by the available observations. To increase the 
signal-to-noise ratio we stacked the last two AES epochs (days 
84 and 86). 

In Table |4] we report the EWs of Nai D2 and Di directly 
measured on all our spectra, together with their associated 
RMS uncertainties, the signal-to-noise ratio achieved at 5900 
A, and the FWHM resolution.The weighted average EWs com- 
puted over all epochs for the two lines are 47.1 ±0.8 mA and 
27.4+1.1 mA, respectively (the standard deviation of the EW 
measurements is 2.6 and 3.4 mA, respectively). The EWs are 
plotted as a function of time in Fig. [5] Although all measure- 
ments are within less than 3-cr from the weighted average (and 
are therefore formally consistent with a null variation), there is 
a systematic increase in both components. While the statistical 
significance for the Dj component is low, the Pearson correla- 
tion factor is 0.82 for D2. A least squares fit to the data gives 
EW(Nai D 2 )=44.6+1.4 + (0.159±0.033) f, where t is the time 
from explosion expressed in days. The corresponding slope for 
D[ is 0.027±0.045 mA day -1 , which is statistically compati- 
ble with a null variation. However, we notice that in conditions 
of mild saturation (as is the case here), EW(D2)/EW(Di)x:1.9. 
Therefore, assuming the measured slope for the D2 component, 
the one expected for Di is «0.08 mA day 1 , which is not in- 
compatible with the low measured value. 

Considering the two extreme values (days 15.1 and 86.4), 
the observed peak-to-peak Nai D2 EW variation amounts 
to 15.6±6.5 mA (considering days 15.1 and 85.2 this is 
14.4±4.1 mA, exceeding the 3-cr level). Although statistically 
significant, this is very small when compared to the variations 
seen in SN 2006X (Patat et al. |2007at , and SN 20071e (Simon 
et al. 2009) during the first 3 months of their evolution, which 
exceeded 100 mA0 



3 Even larger variations were reported for SN 1999cl (Blondin et al. 
[2009b and SN 2006dd (Stritzinger et al. [20TDl . but no high resolution 
spectroscopy was available in those cases. 



6. Implications on ISM small-scale structure 

Given the high velocity of the SN ejecta (>10 4 kms -1 ), the 
photosphere of a Type la SN expands very rapidly, reaching a 
radius larger than 100 AU (~10 15 cm) at maximum light, i.e. 
about 20 days after the explosion. Because of this, the light 
an observer receives from the SN samples different regions of 
intervening inter-stellar material (ISM) as time goes by. If the 
ISM is patchy on the scales of a few tens of AUs, and the den- 
sity fluctuations are large enough, this can translate into vari- 
able absorption features. This geometrical effect and its appli- 
cation to the study of ISM small scale structure are presented 
in Patat et al. (2010), to which we refer the reader for a detailed 
discussion. If we include the EW measurement published by 
Nugent et al. (2011c), the available data allow us to apply this 
method for the first three months of the SN evolution. During 
this time interval the photospheric radius increased from ~10 
AU to ~350 AU (the photospheric radius at our first epoch is 
-60 AU. See Patat et al. EoTOt . 

From the fact that the Nai changed by <15 mA over three 
months, we conclude that the intervening material generating 
the feature at ~180 kms~' cannot be confined in an isolated, 
small (r<200 AU) cloud. As the simulations show (Patat et al. 
120101 their Fig. 3), such small clouds with central column den- 
sities of a few 10 1 1 cirT 2 would be easily detected with the time 
coverage and the EW accuracy characterizing our data set. If 
we adopt an upper limit of 15 mA for the variation of EW(D2), 
we conclude that the cloud responsible for the absorption at 
~180 km s _1 must have a radius larger than 250 AU. This does 
not come as a surprise. The existence of small isolated clouds 
(r~100 AU, N(Nai)~10 n -10 12 cirT 2 ), in the past largely in- 
voked to explain the observed Hi column density fluctuations, 
is severely questioned in terms of pressure equilibrium argu- 
ments, and the yet unknown mechanisms that would produce 
them (see Heiles [19971 and references therein). It is therefore 
more physically meaningful to consider a diffuse ISM, with a 
fractal structure described by a single power law. As it has been 
shown by Deshpande (2000), this scenario provides a convinc- 
ing match to the observations down to AU scales. 

The high-resolution coverage of SN 201 lfe offers a unique 
opportunity to derive the small scale structure properties of 
extra-galactic ISM, a task which can be hardly addressed with 
any other known technique. For this purpose we have com- 
puted a series of models along the lines described in Patat et 
al. (j2pl0 ), simulating the observations of our data set, to which 
we have added the very early epoch by Nugent et al. (1201 lcl) . 
The average column density was set to N(Nai)=2.8xlO n cirT 2 
(corresponding to the weighted average EW(Na 1 D2)), and the 
velocity parameter Z?=3.8 km s _1 , as deduced from the line pro- 
file fittings. For the power-law spectral index we have used the 
canonical value y=-2.8 (Deshpande 2000 ). Given the stochas- 
tic nature of the fractal clouds, we have adopted a statisti- 
cal approach, in which the time dependence of EW(Na 1 D2) 
is examined for a large number (5000) of cloud realizations, 
as a function of the peak-to-peak variation AN on the scale 
of 1000 AU. The simulations show that peak-to-peak varia- 
tions AEW <15 mA imply AN <4xlO n cirT 2 (AN/N £1.4). 
With this density contrast, maximum epoch-to-epoch varia- 
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Fig. 5. EW time evolution of Nai D2 (upper panel) and Di 
(lower panel). Values are from Table [4] The short-dashed hori- 
zontal lines mark the weighted average values, whilst the long- 
dashed lines are least squares fits to the data. 



tions reach 9 mA between the first (day 1.5) and the sec- 
ond epoch (day 8), with an average absolute variation of 
~2 mA. The maximum epoch-to-epoch variations decrease be- 
low ~3 mA already as of day 1 1, as the smoothing effect of the 
extended SN photo-disk grows with time. The maximum devia- 
tion measured in our simulations between any pair of epochs in 
the simulated data-set is ~16 mA, while the 99-th percentile of 
the distribution is 12 mA. The corresponding average absolute 
variation is ~4 mA, with a semi-interquartile range of 2 mA. 
The RMS EW variation from the weighted average is 2.6 mA, 
that matches the value derived from the data (2.7 mA, see pre- 
vious section). The simulated EW(t) curves were analyzed in 
detail, with the aim of understanding whether the observed 
time dependency is consistent with the fractal ISM model. They 
show a smooth behavior, with the EW steadily increasing or de- 
creasing, depending on whether the SN appears projected onto 
an under- or an over-density region. Slopes |dEW/dt| >0. 1 mA 
day -1 are observed in ~1 1% of the cases for AN/N =1.4, while 
the fraction increases to ~24% for AN/N -2. Therefore, the 
relatively large slope derived from the observations is some- 
what disfavored with respect to milder time dependencies, but 
still consistent with the adopted model for the ISM structure. 
For comparison, SN 20071e showed a number of Nai com- 
ponents that remained constant to within a few mA for sev- 
eral months (Simon et al. 120091 ), implying a null variation 
(i.e. |dEW/dt| <0.03 mA day -1 ). Our simulations show that 
for a column density N(Nai)=3xl0 n cirT 2 and AN/N-1.4 
this occurs in ~39% of the realizations. For larger fluctuations 
(AN/N=2) this probability decreases to ~28%. Thus, the two 
different behaviors observed in SNe 201 lfe and 20071e are not 
in contradiction, although the latter is statistically more favored 
than the former. 



The derived density contrast is compatible with the values 
measured in the Milky Way (Frail et al. |1994| l and, therefore, 
the upper limit on the column density variation is fully in line 
with what is expected for a diffuse ISM having the same prop- 
erties as in our own Galaxy. In other words, even in the case the 
observed fluctuations were real (including the apparent steady 
growth), these are consistent with the geometrical effect de- 
scribed in Patat et al. (|2010). Therefore, we conclude that the 
absorption feature observed at ~180 km s _I is most likely asso- 
ciated to the ISM, and is placed at distances larger than ~10 19 
cm (~10 light years) from the explosion site, where the effects 
of the SN radiation field on the Na 1 ionization balance are neg- 
ligible (Patat et al. |2007bl Simon et al. 120091 . Unfortunately, 
the signal-to-noise ratio at ~3900 A in our data is not suffi- 
cient to study with the required accuracy the evolution of the 
Ca 11 H&K features, which in SNe 2006X and 20071e showed a 
distinct behavior (Patat et al. l2007al Simon et al. 12009b . 

Because of the relatively small distance, it will be possible 
to follow SN 201 lfe with high-resolution spectroscopy as the 
object fades away. Although this will require larger telescopes, 
the data may show further variations. As the SN is now entering 
its nebular phase, it is difficult to make predictions based on our 
model, because this is based on the presence of a photosphere 
(Patat et al. I201QI >. However, as this recedes into the ejecta and 
the geometry of the emitting material changes with time, it is 
possible that similar effects take place. 



7. Discussion and conclusions 



In their analysis, Nugent et al. d201 let conclude that SN 201 lfe 
is a slightly sub-luminous, Branch-normal explosion, very sim- 
ilar to the well studied SN 1994D (Patat et al. |1996) . However, 
it is by no means clear whether all "normal" events, sharing 
very similar spectroscopic and photometric properties, do also 
share a common progenitor channel. For this reason it is impor- 
tant to tie the [reasonably well] known progenitor nature of this 
particular event to other properties, which may allow us to iden- 
tify similar events in less favorable cases, i.e. at larger distances 
and with no such exceptionally early discoveries. Among these 
is the circumstellar environment. 

In the last five years, an alternative way of probing the ex- 
plosion environment has been deployed. The method and its 
first application were presented in Patat et al. (2007a), and fur- 
ther expanded in Simon et al. (2009), to which we refer the 
reader for a more quantitative treatment. The idea is based on 
the fact that, at variance with core-collapse events, SNIa are 
relatively weak UV sources. As a consequence, their ability 
to ionize possible circumstellar matter is rather limited in dis- 
tance. Therefore, provided that the gas is placed along the line 
of sight, it can be revealed by optical absorption lines like Ca 11 
H&K, Na 1 D, and K 1. These are all strong lines, and hence they 
enable the detection of tiny amounts of material. As opposed to 
the emission lines that are expected to arise only in the case of 
direct ejecta-CSM interaction (mainly H and He), absorption 
lines are generated by material placed virtually at any distance 
along the line of sight. 



10 



F. Patat et al.: Multi-epoch, high-resolution spectroscopy of SN 201 lfe 



This method has now been applied to a number of Type la 
SNe, and in two cases has led to the detection of time-variant 
Nai features, which were interpreted as arising in the mate- 
rial lost by the progenitor system before the explosion (Patat 
et al. I2007al Simon et al. 12009 ). The velocities of this mate- 
rial are consistent with those of a red giant, hence favoring a 
symbiotic system. The complex velocity structure shown by 
SN 2006X was explained as being the possible outcome of 
recurrent-novae episodes (Patat et al. 2007a). This suggestion 
is supported by the results of a similar analysis carried out on 
the known recurrent nova RS Oph, which shows strikingly sim- 
ilar CSM properties (Patat et al. 1201 11) . A study conducted on 
a large number of SNIa has shown a statistically significant 
excess of blue-shifted Nai features, which were interpreted as 
arising in material lost by the progenitor system (Sternberg et 
al. 1201 11 1. This has given support to the conclusions reached for 
SNe 2006X and 20071e, which have been suggested to be con- 
nected to binary systems hosting a mass-losing giant. If this in- 
terpretation is correct, the measured fraction of SNIa display- 
ing blue-shifted features suggests this channel is not a negligi- 
ble contributor to the global rate of thermonuclear events. This 
conclusion is apparently inconsistent with the limits set by the 
lack of ejecta-wind signatures, and the failed detection of inter- 
action between the ejecta and the companion atmosphere dur- 
ing the very early epochs. Which may simply mean that the 
picture is not complete yet; see for instance the discussion in 
Chomiuk et al. (1201 II ) about our understanding of the wind ge- 
ometry in SNIa progenitors (see also Patat l20lTT ). 

What makes SN 201 lfe special is that the progenitor sys- 
tem has been identified with a significant confidence through 
a number of very stringent constraints (Li et al. 1201 11 Brown 
et al. I20TT1 Nugent et al. 1201 let Bloom et al. I20TT1 Liu et al. 
1201 11 Horesh et al. 1201 11 ). These constraints clearly rule out 
symbiotic systems with a red giant companion, hence exclud- 
ing recurrent novae like RS Oph. At least for this one event. 

The absence of time-variant, blue-shifted features and the 
"clean" environment reported in this paper is fully in line with 
the progenitor nature derived from the very early observations. 
The bulk of the material transferred from the companion re- 
mains confined to within the binary system, and there are no 
recurrent outbursts that would create the complex CSM struc- 
ture seen in RS Oph (Patat et al. 1201 II ). and possibly present 
in SN 2006X (Patat et al. I2007al ) and in a significant fraction 
of other SN la (>20%; Sternberg et al. 1201 lb . The key question 
is whether the argument can be reversed, and the absence of 
such features in any SN la can be interpreted in terms of a sim- 
ilar progenitor, at least in a statistical sense. The answer to this 
question is awaiting detailed hydrodynamical modeling, which 
would tell us how frequently the physical conditions for pro- 
ducing observable narrow features are met, and what the view- 
ing angle effects are (Mohamed et al., in preparation). Before 
this is accomplished, we can only speculate based on very small 
numbers. 

As pointed out by Nugent et al. d201 lcb . SN 201 lfe is spec- 
troscopically very similar to the somewhat sub-luminous events 
SN 1992A (Kirshner et al. [19931 ) and SN 1994D (Patat et al. 
1996 ). Although the link between these objects needs to be 
assessed considering a number of observables, here we can- 



not refrain from noticing that all three objects present a very 
"clean" circum/inter-stellar environment. To illustrate this fact, 
in Fig. [6] we compare the high-resolution spectra of SN 1994D 
and SN 201 lfe in the velocity ranges corresponding to the host 
galaxies recession (no measurable absorption features were 
reported for SN 1992A; Kirshner et al. [19931 . Although for 
SN 1994D there is no time coverage (see also King et al. 1995 ), 
the two spectra are remarkably similar. 

In the schema devised by Sternberg et al. (1201 II) . SNe 
1994D and 201 lfe would be classified in the single/symmetric 
group (SS-type, 14.3%), while SN 1992A would fall in the 
group with no absorption detection (N-type, 37.1%). An in- 
spection of the SS group shows that 3 SS-type events (out of 5) 
have a single component, similar to SN 201 lfe. If we assume 
that in all these cases the absorption is attributable to the host 
galaxy (and not to the CSM environment), and we add them 
to the N-type events, then the fraction of "clean events" in the 
Sternberg et al. sample is ~46%. One can imagine cases where 
the CSM environment is "clean", but the line of sight inter- 
cepts ISM clouds that are totally unrelated to the progenitor. 
Therefore, while the lack of conspicuous absorption features 
can be confidently interpreted in terms of a gas -free environ- 
ment, the reverse argument does not hold. 

The projected position of SN 201 lfe falls close to a periph- 
eral spiral arm (see for instance Nugent et al. 1201 let . However, 
the very low reddening and the absence of complex absorp- 
tion structures strongly suggest the explosion took place on the 
front side of the face-on host, most probably well above its 
galactic plane. Had the SN been placed on the rear side of the 
galaxy, the IS imprints would have been much more complex, 
irrespective of the progenitor's nature deduced from the early 
data. Therefore the above fraction can only represent a lower 
limit. On the other hand, because of the viewing angle effects, 
the fraction of events allegedly produced by systems hosting a 
giant estimated by Sternberg et al. ( 1201 II ) is also a lower limit 
(20-25%). 

The sample of SN la with sufficient spectral resolution and 
time coverage is still relatively small, and therefore drawing 
conclusions may be premature. The fraction of events show- 
ing time-variant features is ~17% (Patat [20TTT ). Probing only 
the line of sight, this indicator also suffers from viewing angle 
effects. In addition, the method only reveals material placed 
in the suitable range of distances from the SN, and having 
the required densities (Patat et al. I2007al Simon et al. 120091) . 
Therefore, again, it gives only a lower limit to the fraction of 
events generated through this channel. Interestingly, this limit 
is reasonably similar to that derived from the excess of blue- 
shifted features (20%; Sternberg et al. 120111 ). and it confirms 
the conclusion that this channel contributes in a significant way 
to the global SN la rate. This result has still to come to terms 
with the opposite evidences given by the lack of CSM signa- 
tures in radio, UV, and X-rays domains, and the failed detection 
of extended companion atmospheres in the early stages of the 
SN evolution. 

Although expensive in terms of telescope time, multi- 
epoch, high-resolution spectroscopy of SN la constitutes a 
powerful and independent tool for investigating the nature of 
SN la progenitors. With the future availability of statistically 
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Fig. 6. Comparison between the Nai D spectral regions of 
SN 1994D (upper; Patat et al. [T996l i. and SN 201 lfe (lower; 
SARG). Heliocentric velocities are computed with respect to 
the D2 component. The resolution of the SN 1994D spectrum 
is ~10 kms -1 . 

significant samples, it will most likely help us in settling the 
matter about the contribution of the various channels leading 
to Type la explosions. And all the pieces that currently appear 
to form a contradictory set of information may fall together to 
form a consistent picture. 
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